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ABSTRACT 

We examine the effect of Lyman continuum photons on the 21 cm background in the high- 
redshift universe. The brightness temperature of this transition is determined by the spin tem- 
perature T s , which describes the relative populations of the singlet and triplet hyperfine states. 
Once the first luminous sources appear, T s is set by the Wouthuysen-Field effect, in which 
Lyman-series photons mix the hyperfine levels. Here we consider coupling through n > 2 
Lyman photons. We first show that coupling (and heating) from scattering of Lyri photons is 
negligible, because they rapidly cascade to lower-energy photons. These cascades can result 
in either a Lya photon - which will then affect T s according to the usual Wouthuysen-Field 
mechanism - or photons from the 2s — > Is continuum, which escape without scattering. We 
show that a proper treatment of the cascades delays the onset of strong Wouthuysen-Field 
coupling and affects the power spectrum of brightness fluctuations when the overall coupling 
is still relatively weak (i.e., around the time of the first stars). Cascades damp fluctuations on 
small scales because only ~ 1/3 of Lyn photons cascade through Lya, but they do not affect 
the large-scale power because that arises from those photons that redshift directly into the Lya 
transition. We also comment on the utility of Lyn transitions in providing "standard rulers" 
with which to study the high-redshift unvierse. 
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1 INTRODUCTION 

One potentially promising probe of the co smic dark ages is 21 
cm tomography. It has long been known faogan&Ree3ll979t 
IScott & Reesl Il99d) that neutral hydrogen in the intergalactic 
medium (IGM) may be detectable in emission or absorption against 
the cosmic microwave background (CMB) at the wavelength of 
the redshifted 21 cm line, the spin-flip transition between the sin- 
glet and triplet hyperfine levels of the hydrogen ground state. 
The brightness of this transition wil l thus trace the dis tribution 
of HI in the high-redshift universe (Field 1958, 1959a), which 
gives the signal angular structure as well as structure in redshift 
space. These features arise from inhomogeneities in the gas density 
field, the hydrogen ionization fraction and the spin temperature. 
Mad au. Meiksin & Reesl (1997) showed that the first stars could 
cause a rapid evolution in the signal through their effect on the spin 
temperature. Consequently, the 21 cm signal can provide unparal- 
leled information about the "twilight zone" when the first luminous 
sources formed and the epoch of reionization and reheating com- 
menced. 

Despite the theoretical promise of this probe, it is only 
with improvements in computing power that building radio ar- 
rays with sufficient sensitivity, capable of correlating billions of 
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visibi lity measurements, has become possible jMorales & Hewittl 
12004 . Three such arrays (LOFAR 1 , MWA 2 , and PAST 3 ) will 
soon be operational, opening a window onto this new low fre- 
quency band. Before a detection can be made, however, there are 
still major scientific and technical challenges to be met. Iono- 
spheric scattering and terrestrial interference are two serious is- 
sues. Also worrying is the need to remove foregrounds, which are 
many orders of magnitude stronger than the signal. Multifrequenc 
subtraction techniques ( Zaldarriaaa. Furlanetto & Hernauist 200- 

sir 



icy 
)4; 



iMorales & Hewitt 2004 ISantos. Coorav ~ Knox 2005), exploiting 
the smoothness of the foreground spectra, have been proposed, but 
their effectiveness has yet to be tested. The challenges are great, but 
so are the opportunities. It is thus crucial to understand the nature 
of the 21 cm signal as we commence these searches. 

Fluctuations in the 21 cm signal arise from both cosmo- 
logical and astrophysical sources. Most previous work has fo- 
cusse d on the signal due to den sity perturbations I Mad au et alJ 
1 19971 |Loeb_&_Zaldjirriag3 |200 4) or from inhomogeneous ion- 
ization j Oardi & Madatj]2003|jFurlanetto, Sokasian & Hernauist 
2004 ; Furlanetto. Zaldarriaga & Hernauist 2004J). An additional 



1 See http://www.lofar.org/ 

2 See http://web.haystack.mit.edu/arrays/MWA/ 



See Pen et al. (2005). 
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source of fluctuations is the spin temperature, which describes the 
relative occupation of the singlet and triplet hyperfine levels. These 
levels may be excited by three primary mechanisms: absorption of 
CMB photons, atomic collisions, and absorptio n and re-em ission 
of Lya phot ons (the Wouthuysen-Field effect; IWouthuvserlll952t 
Field 1959a). The first two processes rely upon simple physics, but 
the last one allows us to study the properties of luminous sources, 
which determine the background radiation field. 



iBarkana & Loebl l2005bl henceforth BL05) studied the signal 
generated by the first generation of collapsed objects. These high 
redshift objects are highly biased, leading to large variations in their 
number density. This, combined with the 1 /r 2 dependence of the 
flux, causes large fluctuations in the Lya background, which can be 
probed through their effect on the 21 cm transition. Exploiting the 
anisotropy induced by peculiar velocities (Bharadwai & Ali 2004; 
IBarkana & Loeb 2005a), they showed that information about the 
Lya radiation field could be extracted from the power spectrum 
of 21 cm fluctuations and separated into those fluctuations corre- 
lated and uncorrelated with the density field. The features of these 
spectra allow extraction of astrophysical parameters such as the star 
formation rate and bias. However, it is not a trivial task to relate 
the emissivity to a distribution of Lya photons. The background in 
the Lya line is composed of two parts: those photons that have red- 
shifted directly to the Lya frequency and those produced by atomic 
cascades from higher Lyman series photons. To calculate this latter 
component, BL05 assumed that atomic cascades were 100% effi- 
cient at converting photons absorbed at a Lyman resonance into a 
Lya photon, while in reality most cascades end in two photon de- 
cay from the 2S level. 

In this paper, we calculate the exact cascade conversion prob- 
abilities from basic atomic physics. In addition, we discuss the pos- 
sibility of level mixing by scattering of Lyn photons via a straight- 
forward generalisation of the Wouthuysen-Field effect. We then ap- 
ply the cascade efficiencies to calculate the Lya flux profile of an 
isolated source. The existence of discrete horizons, determined by 
the maximum distance a photon can travel before it redshifts into a 
given Lyman resonance, imprints a series of discontinuities into the 
profile, which can in principle be used as a standard ruler. We apply 
these results to the power spectra of 21 cm fluctuations during the 
epoch of the first stars, showing that these corrections can not be 
ignored when extracting astrophysical parameters. 

The layout of this paper is as follows. In S|2|we introduce the 
formalism for describing 2 1 cm fluctuations and the dominant cou- 
pling mechanism, the Wouthuysen-Field effect. In J3]we discuss 
the possibility of direct pumping by Lyn photons. Next, in Sj4] 
we detail the atomic physics of radiative cascades in atomic hy- 
drogen. The results are applied to the Lya flux profile of an iso- 
lated source in fj5]and to the 21 cm power spectrum from the first 
galaxies in We also discuss some of the limitations of this for- 
malism. Finally, we summarise our results in In an Appendix, 
we review the equations needed to calculate analytically the Ein- 
stein A coefficients for the hydrogen atom. Throughout, we assume 
(fi m , fii,, n A , h, as, n s ) = (0.3, 0.046, 0.7, 0.7, 0.9, 1.0), c onsis- 
tent with the most recent measurements ISpergel et alj2003l) . 

During the preparation of this paper, iHiratal 120051) sub- 
mitted a preprint covering similar material. We have confirmed 
agreement where there is overlap. The main results of this 
work were discussed at the "Reionizing the Universe" confer- 
ence in Groningen, The Netherlands (June 27- July 1, 2005; see 
http://www.astro.rug.nl/~cosmo05/program.html I. 



2 21 CM FORMALISM AND THE 

WOUTHUYSEN-FIELD MECHANISM 

The 21 cm line of the hydrogen atom results from hyperfine split- 
ting of the IS* ground state due to the interaction of the magnetic 
moments of the proton and the electron. The HI spin temperature 
T s is defined via the relative number density of hydrogen atoms in 
the 15 singlet and triplet levels ni/no = (31/30) exp(— T* /T s ), 
where (g± /go) — 3 is the ratio of the spin degeneracy factors of the 
two levels, and T* = /ic/fcA2icm = 0.0628 K. The optical depth 
of this transition is small at all relevant redshifts, so the brightness 
temperature of the CMB is 

where the optical depth for resonant 21 cm absorption is 

3c\ 2 hAionm 
T ~ 327rfcgT s (l + z)(dv r /dr 

Here nm is the number density of neutral hydrogen, A\$ — 2.85 x 
10 -15 s _1 is the spontaneous emission coefficient, and dv r /dr is 
the gradient of the physical velocity along the line of sight with r 
the comoving distance. When T 3 < Tomb there is a net absorption 
of CMB photons, and we observe a decrement in the brightness 
temperature. 

The spin temperature is determined by three coupling mecha- 
nisms. Radiative transitions due to absorption of CMB photons (as 
well as stimulated emission) tend to drive T s — * Tomb- Spin flips 
from atomic collisions drive T s — + Tk, the gas kinetic tem pera- 
ture. F inally, the Wouthuysen-Field effect 1 Wout huvsenll952HFieldl 
1958), which is the main focus of this paper, also drives T s — > Tk 
(see below). The combination that appears in Q can be written 

T s — Tomb 



(2) 



Xtot 



1 - 



Tomb 1 



(3) 



T a 1 + xtot V Tk 

, + x c is the sum of the radiative and collisional 



where x to t = x 
coupling parameters. The latter is 

_ 4 K i_o(T fe )n H T, 



SAioTcmb ' 

where ki_q is tabulate d as a function of Tk lAllison & Da lgarno 
1 19691 IZvgelmarJl2005l) . The spin temperature becomes strongly 
coupled to the gas temperature when x to t > 1. 

A schematic diagram of the Wouthuysen-Field effect is shown 
in FigureQ it mixes the hyperfine levels through absorption and re- 
emission of Lya photons. Quantum selection rules allow transitions 
for which the total spin angular momentum F changes by AF = 
0, ±1 (except — > 0), making only two of the four n = 2 levels 
accessible to both the n = 1 singlet and triplet states. Transitions 
to either of these states can change T a . The coupling coefficient is 

4P a T* 

Xc = — -. , (5) 

27AioTcmb 

where P a is the Lya scattering rate !Mada uetalll997l) . If resonant 
scattering of Lya photons occurs rapidly enough T 3 will be driven 
to T a , the colour temperature of the radiation field at the Lya fre- 
quency jFieldll958t iMadau et alJll997l) . In parallel, the repeated 
scattering of Lya photons by t he thermal distribution of atoms 
brings T a -> T k iFieldHl959bl lHiratdl2005l) . Consequently, the 
Wouthuysen-Field effect provides an effective coupling between 
the spin temperature and the gas kinetic temperature. 

We can also write the Wouthuysen-Field coupling as 

167r 2 T*e 2 /, 



27 AioTcMsm e c 



Set J a 



(6) 
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Figure 1. Hyperfine structure of the 2P and 15 level of the hydrogen atom. 
Levels are labelled according to the notation n pLj, where n, L, and J are 
the usual radial, orbital angular momentum and total angular momentum 
quantum numbers. F = I + J is the quantum number obtained from the 
nuclear spin and J. Allowed transitions obey AF = 0, ±1 (except — » 
0). Those relevant for the Wouthuysen-Field effect are indicated by solid 
curves, while dashed curves indicate the remaining allowed transitions. 



where f a — 0.4162 is the oscillator strengt h for the Lya tran sition, 
S a is a correction factor of order unity l Chen & Miralda-Escude 
|2004 lHiraill2005h that accounts for the redistribution of pho- 
ton energies due to repeated scattering off the thermal distribu- 
tion of atoms, and J a is the angle-averaged specific intensity of 
Lya photons by photon number. For reference, a Lya flux of 
J a = 1.165 x 1(T 10 [(l + z)/20]cm- 2 s -1 Hz -1 sr -1 yields 
x a = 1 (corresponding to P a = 7.85 x 1CT 13 [1 + z] s" 1 ). 

Fluctuations in the brightness temperature arise from fluctu- 
ations in the density, the Wouthuysen-Field coupling, the neutral 
fraction xm and the radial velocity component. To linear order 



; 05 + -rr^-5 Xa + 5 XH1 — 5d r v r , 
Xtot 



(7) 



where 5 a is the fractional perturbation in a, 5 is the fractional 
density perturbation, and x to t ~ Xtot(l + %tot)- P is a param- 
eter describing the thermal history of the gas, whic h we assume 
to hav e cooled adiabatically, so that P w 0.2. iNaoz & Barkanl 
(2005) showed that j3 slightly increases and exhibits mild scale de- 
pendence, as gas temperature fluctuations do not exactly track the 
density fluctuations. We note that, on scales 0.01 Mpc -1 < k < 
10 3 Mpc -1 , P is approximately constant at z = 20 and choose 
to ignore this subtlety for ease of comparison with BL05. The first 
three components of equation Q are isotropic, but the velocity fluc- 
tuation introdu ces an anisotropy of the form 5d T v T (k) = — u 2 <5 in 
Fourier space iBharadwai & Alill2004l) . where u is the cosine of 
the angle between the wavenumber k of the Fourier mode and the 
line of sight. This allows us to separate the brightness temp eratu re 
power spectrum Pr b into powers of u 2 iBarkana & Loeb 2005a) 



P Tb (k) = M 4 P M 4 (fc) + tfPp (fc) + P^O (k). 



(8) 



depends only on the matter power spectrum. P^i contains cross- 
correlations between matter fluctuations and both 5 X a and 5 XUI , 
making it an ideal probe of fluctuations in the radiation background. 
In particular, at sufficiently high redshifts such that xhi <C 1, it 
probes variations in the Lya background. Linear combinations of 
these three terms can be used to extract detailed information about 
other types of fluctuations iBarkana & Loebl2005al) . 



3 DIRECT PUMPING BY LYMAN SERIES PHOTONS 

Of course, the radiation background contains photons that redshift 
into all the Lyman transitions, not just Lya. The main purpose of 
this paper is to examine how these affect T a . The existing litera- 
ture assumes that all Lyra photons are immediately converted into 
Lya photons by atomic cascades (e.g. BL05). In reality, there are 
two different contributions to consider: one due to scattering of the 
Lyra photon itself and the other due to its cascade products. In this 
section, we discuss the direct contribution of Lyra scattering to the 
coupling of T s and Tk , which occurs in a manner exactly analogous 
to the Wouthuysen-Field effect. For this effect to be significant two 
requirements must be fulfilled. First, the scattering rate of Lyra pho- 
tons must be sufficient to couple T 3 and T n , the Lyra colour tem- 
perature. Second, it must be sufficient to drive T n — > Tk. We will 
argue that neither condition is satisfied in practice. 

The IGM is optically thick r > 1 to all Lyman series tran- 
sitions with ra < 100. Consequently, a L ya photon emitted by a 
star will scatter many times (~ r ~ 10 6 : lGunn & Peterson! 19651) 
before it finally escapes by redshifting across the line width; each 
of these scatterings contributes to the Wouthuysen-Field coupling. 
A Lyra photon can escape by redshifting across the line width, but 
a transition to a level other than ra = 1 will also remove it. The 
probability for a decay from an initial state i to a final state / is 
given in terms of the Einstein Ai f coefficients by 



Pi/ = 



Ai 



f 



(9) 



Appendix [A] summarises the expressions needed to compute 
the Einstein Aif coefficients. For the Lyman series transitions 
P n p^is w 0.8 (see Table 1) so that a Lyra photon will scatter 
of order -/V sca t ~ 1/(1 — P n p^is) ~ 5 times before undergoing 
a cascade. 

Because a cascade occurs long before escape via redshifting, 
the coupling from direct pumping is negligible. Recall that the scat- 
tering rate P x for the ph oton type X =Lya, Ly/3, etc. may be 
expressed as iFiel dll959al) 



nmPx 



cat TlX 



(10) 



The anisotropy is sourced only by density fluctuations, so that P M 



in terms of the production rate of photons per unit volume hx- 
It is then clear that, for similar production rates (i.e., for sources 
with a reasonably flat spectrum), P n /P a ~ iV sca t jn /iV sca t,a ~ 
5 x 10~ 6 . This simple argument shows that the contribution from 
direct pumping by Lyra photons will be negligible compared to that 
of the Lya photons, because x a oc P a . 

The second question, whether T n — » Tk, is still relevant 
for heating of the gas by repeated scatterings. Given the reduced 
number of scattering events, it seems unlikely to be the case, 
but a full c alculation using a Mo nte Carlo method or follow- 
ing C hen & Miralda-Escudel 120041) is required to rigourously an- 
swer this question. Lack of equilibrium would make Lyra scat- 
tering a more efficient source of heat, on a per scattering basis. 
IChen & Miralda-Escudel j2004) have shown that Lya heating is 
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much smaller than previous calculations indicated I Mad au et alJ 
Il997t) . because T a ~ Tk, which reduces the heat transferred per 
collision. This is unlikely to be the case for the Lyn. 

Following Madau et al ill 9971) . we can estimate the maximum 
heating from a single Lyn scattering by assuming that all of the 
atomic recoil energy for a stationary atom is deposited in the gas. 
Momentum conservation then demands 



E n 



AE 

E 



hv n P n , 



(11) 



where (AE / E) ~ 10~ 8 is the fraction of energy lost by a Lyn 
photon after scattering from a stationary hydrogen atom, and hv n is 
the energy of the photon. Assuming the production rate of Lyn pho- 
tons is comparable to that of the Lya photons and taking x a — 1, 
we then obtain E n ~ 0.002[(1 + z)/10] KGyr -1 . This is much 
smaller than the Lya heating rate, even including the T a ~ Tk cor- 
rection, so we do not expect Lyn scattering to be a significant heat 
source. Furthermore, if T n ~ Tk the rate would be mu ch smaller 
than this estimate, as in IChen & Miralda-Escudd <2004 . 




n=l 



D 



4 LYMAN SERIES CASCADES 

An excited state of hydrogen may reach the ground state in three 
ways. Firstly, it may decay directly to the ground state from an nP 
state (n > 2), generating a Lyn photon. Secondly, it may cascade 
to the metastable 2S level. Decay from the 25* level proceeds via 
a forbidden two photon process. Finally, it may cascade to the 2P 
level, from which it will produce a Lya photon. We are primar- 
ily interested in the fraction of decays that generate Lya photons, 
which will increase the Lya flux pumping the hyperfine levels. 

The fraction of cascades that generate Lya photons can be de- 
termined straightforwardly from the selection rules and the decay 
rates. As an example, consider the Ly/3 system. Absorption of a 
Ly/3 photon excites the atom into the 3P level. As illustrated in 
Figure [2] the 3P level can decay directly to the ground state, re- 
generating the Ly/3 photon, or to the 2S level, where it will decay 
by two photon emission. The selection rules forbid Ly/3 photons 
from being converted into Lya photons. In contrast, the 4P level, 
excited by absorption of Ly7, can cascade via the 35* or 3D levels 
to the 2P level and then generate Lya. 

To calculate the probability /recycle that a Lyn photon will 
generate a Lya photon, we apply an iterative algorithm. The ex- 
pression 

/ 'recycle,/ — ^ ^ Pi/ /recycle,/ 

(12) 

/ 

relates the conversion probability for the initial level i to the con- 
version probabilities of all possible lower levels /. The decay prob- 
abilities are calculated using equation {5). We then iterate from low 
to high n, calculating each /recycle in turn. 

In our particular case of an optically thick medium, we can 
ignore direct transitions to the ground state. These generate a Lyn 
photon, which will rapidly be reabsorbed and regenerate the nP 
state. Therefore, such decays will not affect the net population of 
photons or of excited states. We incorporate this into the calculation 
by setting A n p^\s ~ I Furlane tto et all2005l) . 

Results for the lowest Lyman series transitions are sum- 
marised in Table 1 and plotted in Figure |5] 4 These results are in 

4 Code for calculating the conversion factors is available at 
http://www.tapir.caltech.edu/~jp/cascade/ 



Figure 2. Energy level diagram for the hydrogen atom illustrating Ly/3 and 
Ly7 cascades. Marked decays are distinguished as cascades (solid curves), 
Lyn (dashed curves), Lya (dot-dashed curve), and the two photon decay 
(dotted curve). Note that the selection rules (AL = ±1) decouple the 3P 
and 2P levels, preventing Ly/3 from being converted into Lya. 
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0.8817 


18 


0.3561 


0.7748 
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0.2609 


0.8390 


19 


0.3565 


0.7743 
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0.3078 


0.8178 


20 


0.3569 


0.7738 


6 


0.3259 


0.8053 


21 


0.3572 


0.7734 


7 


0.3353 


0.7972 


22 


0.3575 


0.7731 


8 


0.3410 


0.7917 


23 


0.3578 


0.7728 


9 


0.3448 


0.7877 


24 


0.3580 


0.7725 


10 


0.3476 


0.7847 


25 


0.3582 


0.7722 


11 


0.3496 


0.7824 


26 


0.3584 


0.7720 


12 


0.3512 


0.7806 


27 


0.3586 


0.7718 


13 


0.3524 


0.7791 


28 


0.3587 


0.7716 


14 


0.3535 


0.7780 


29 


0.3589 


0.7715 


15 


0.3543 


0.7770 


30 


0.3590 


0.7713 



Table 1. Recycling fractions /recycle an d decay probabilities to the ground 
state, P nP ^ ls . 

agreement with those of Hirata 1 2005). At large n, the conversion 
fractions asymptote to /recycle ~ 0.36 because nearly all cascades 
pass through lower levels. We emphasise again that the quantum 
selection rules forbid a Ly/3 photon from producing a Lya photon. 

Finally we briefly comment on the two p hoton decay f rom the 
25* level (for a more detailed discussion, see Hirata 120051) . Selec- 
tion rules forbid electric dipole transitions from the 2S level to the 
ground state, but the second order two photon decay process can 
occur with A 77 = 8.2 s~ x < A 2 p->is- At z < 400 the CMB 
flux density is sufficiently small that radiative excitations from the 
25* level are negligible. Additionally, at the relevant densities col- 
lisional excitatio n to the 2P level i s slow compared to the two 
photon process (Breit & Teller 1940). Consequently, the 25 level 
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Figure 3. Recycling fractions for Lyn photons. Note that the values level 
off at / rcC y C i c ~ 0.359 and that none of the photons incident on the Ly/3 
resonance are converted into Lyct photons. 



will preferentially decay via this two photon process. These transi- 
tions may themselves affect T s , because both the 2S and IS levels 
have hyperfine structure and any imbalance in the decay constants 
here would affect the IS populations. However, even without de- 
tailed calculations, we can see that the resultant coupling must be 
small. Cascades that do not generate Lya must reach the 2S level, 
so the fraction of Lyn photons that undergo two photon decay is 
/ 77 (n) = 1 — /rccycie(n) ~ 0.64 (See Table 1). Each such decay 
has A^cat.77 = 1 because the resulting photons are not reabsorbed. 
This is much smaller than jV SC at,c< ~ 10 6 . Consequently, only if the 
coupling per scattering were many orders of magnitude larger for 
two photon decay than for Lya scattering could this effect be sig- 
nificant. 



5 THE Lya COUPLING AROUND A SOURCE 

We can see the effects of these recycling fractions on the Lya cou- 
pling by considering the life of a photon emitted from a given 
source. The photon initially propagates freely, redshifting until it 
enters a Lyn resonance. Because the IGM is so optically thick, the 
photon will then scatter several times until a cascade converts it 
into a Lya photon or two 25 — > IS photons. In the latter case, the 
photons escape to infinity; in the former case, it scatters ~ r times 
before redshifting out of the Lya resonance. This establishes a se- 
ries of closely-spaced horizons, because a photon entering the Lyn 
resonance at z must have been emitted below a redshift 



1 + Zmax(n) = (1 + z) 



[l-(n + l)- 2 ] 
(1-n- 2 ) ' 



(13) 



The number of Lyn transitions contributing Lya photons is thus 
a function of the distance from the source. These horizons imprint 
well-defined atomic physics onto the coupling strength by introduc- 
ing a series of discontinuities into the Lya flux profile of a source. 

Thus the Lya flux, J a , arises from a sum over the Lyn lev- 
els, with the maximum n determined by the distance. The sum is 



ultimately truncated at n max ~ 23 to exclude levels for which the 
horizon lies within the HII region of a typical (isolated) galaxy, as 
only neutral hydrogen contributes to 21 cm absorption (BL05). The 
average Lya background is thus 



= E 



dz ^/recycle (^) 



(1 + -) 2 



4tt H 



(14) 



where i>' n is the emission frequency at z corresponding to absorp- 
tion by the level n at z 



'•(! + *)' 



(15) 



and e(v, z) is the comoving photon emissivity (defined as the num- 
ber of photons emitted per unit comoving volume, per proper time 
and frequency, at frequency v and redshift z). To calculate e(v, z), 
we follow the model of BL05 



e(v,z) = n°f*-^;F sal (z)e b (i>), 



(16) 



96 90, and 6520 respective ly (BL05). In calculating F ga i, we use 
fhe lSheth & Tormenl ll999h mass function dn/d m, which matches 



where n£ is the cosmic mean baryon number density today, /, is 
the efficiency with which gas is converted into stars in galactic ha- 
los (and with which Lyman-continuum photons escape their hosts), 
and F ga ,i(z) is the fraction of gas inside galaxies at z. We model 
the spectral distribution function of the sources tb{v) as a separate 
power law tbiy) oc u as _ 1 between Lya and Ly/3 and between Ly/3 
and the Lyman limit. In our calculations, we will assume Popula- 
tion III stars with spectral index a s = 1.29 between Lya and Ly/3, 
normalised to produce 4800 photons per baryon between Lya and 
the Lyman limit, of which 2670 photons are emitted between Lya 
and Ly/3. In contrast, for Population II stars the numbers are 0.14, 
;sp ective l> 
enNl999l) i 

simulations better than the I Press & SchechteJ mass func- 

tion, at least at low redshifts. We assume that atomic hydrogen 
cooling to a viral temperature T v i r ~ 10 4 K sets the minimum halo 
mass. We will normalise /, so that x a = 1 at z — 20; this yields 
/» = 0.16% when we include the correct /recycle- 

To see what fraction of photons from a given source are 
converted into Lya, we integrate et(v) with the proper weight- 
ing by /recycle- We find that /recycle = 0.63, 0.72, and 0.69 for 
a s (Lya — Ly/3) = 1.29, 0.14, and —1.0 respectively (roughly 
corresp onding to Pop. Ill stars, low metallicity Pop II stars, and 
quasars; Zheng et al. 1997). The total flux is significantly less than 
if /recycle = 1, as has been generally assumed before. Thus Lya 
coupling will take place later if the proper atomic physics are in- 
cluded (typically Az > 1 for fixed source parameters; see also 
lHiratal2005l) . Of course, this is only the average value, and around 
a given source there will be a distance dependence. A gas element 
that can only be reached by photons redshifted from below the Ly/3 
resonance will see an effective /recycle = 1. In contrast, a gas ele- 
ment very close to the source will have /recycle ~ 0.36. This will 
be reflected in the brightness temperature power spectrum (see S]6j. 

The Lya flux profile of a galaxy with M ga i = 3 x 10 10 M 
and our fiducial parameters at z = 20 is plotted in Figure [4] 
In our approximation, J a oc M ga i. Thus obtaining x a > 1 at 
r = 10 Mpc requires a galaxy mass of M ga i = 4.2 x 10 12 Mq, 
corresponding to a 14cr fluctuation in the density field. Obviously, 
individual sources do not induce strong Lyman coupling on large 
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Figure 4. Flux profile of a galaxy of mass M = 3 X 10 10 Mq at z = 20 as 
a function of comoving distance r. For comparison, we have plotted the flux 
profile assuming / rccyc i c = 1 (dashed curve), proper atomic physics (solid 
curve) and including only photons with u a < v < ua (dotted curve). Ver- 
tical lines along the lower axis indicate the horizons for the Lyn resonances. 
The horizontal dashed line shows the value of x c at z = 20, illustrating the 
regime where collisional coupling dominates. 



scales. The conversion of photons from Lyn to Lya steepens the 
flux profile beyond the simple 1/r 2 form. Notice that we have nor- 
malised /* for each curve separately, so that x a = 1 at z = 20. 
Because setting / rccyc i c = 1 weights large n transitions more heav- 
ily (and hence small scales), that curve lies below the others at large 
r. The discontinuities occur at the Lyn horizons. In theory, their po- 
sitions yield standard rulers determined by simple atomic physics. 
In practice, the weakness of the discontinuities, and the overlap- 
ping contributions of other nearby sources, makes it unlikely that 
these discontinuities will be observable for an isolated source (see 
also ij6}. Finally, we note that sharp discontinuities only occur if a 
photon undergoes a cascade immediately after entering a Lyn res- 
onance and if the resulting Lya photons redshift out of the Lya 
resonance immediately. The form er is certainly tru e, but the latter 
will affect the shape significantly <Loeb&Rvbickill999l) . 

Additionally, photon horizons affect the radiation heating of 
the gas around a source. The cascades, which result from scatter- 
ing of Lyn photons, deposit some of their radiative energy into the 
kinetic energy of the gas, and their total heating rate differs from 
"continuum " Lya photons because they are injected as line pho- 
tons IChen & Miralda-Escudel [2004 ) . Thus the Lya heating pro- 
file will differ from 1/r 2 . However, Lya heating is typically much 
small er than other sources, so this is unlikely to be important 
IChen & Miralda -Escudel2004l) . 



6 BRIGHTNESS FLUCTUATIONS FROM THE FIRST 
GALAXIES 

In the previous section, we saw that the proper / rccyc i c affects the 
spatial distribution of x a around each source. The most important 
manifestation of this will occur when the Wouthuysen-Field effect 



is just becoming important, around the time of the first galaxies 
(BL05). Those authors showed that Lyn transitions enhance the 
small-scale fluctuations in but they assumed that / rccyc i c = 1. 
We will show how the scale-dependent / rC c y cic modify this signal. 
It is possible to exploit the separation of powers to probe sepa- 
rately fluctuations that correlate with the density field and those, 
like Poisson fluctuations, that do not. We consider each in turn and 
compare to the results of BL05. We set 5 xml = throughout. We 
also assume that the IGM cools adiabatically, with no heat input 
from X-rays. Note that for ease of comparison wit h BL05, we do 
not incorporate the low-temperature corrections of iHiratal 120051) 
(and in any case they are small in our example). 



6.1 Density fluctuations 

Density perturbations source x a fluctuations via three effects 
(BL05). First, the number of galaxies traces, but is biased with 
respect to, the underlying density field. As a result an overdense 
region will contain a factor [1 + b(z)5\ more sources, where b(z) 
is the (mass-averaged) bias, and will have a larger x a . Next, pho- 
ton trajectories near an overdense region are modified by gravita- 
tional lensing, increasing the effective area by a factor (1 + 25/3). 
Finally, peculiar velocities associated with gas flowing into over- 
dense regions establish an anisotropic redshift distortion, which 
modifies the width of the region corresponding to a given observed 
frequency. These three effects may be represented using a linear 
transfer function W(k) relating fluctuations in the coupling S Xa to 
the overdensity S 



5 Xa = W(k)5. 



(17) 



We compute W(k) for a gas element by adding the coupling due 
to Lya flux from each of the Lyn resonances (BL05) 



W(k) = - £ 



,dxi n) D(z') 

dz 



dz' D(z) 

x hi + b(z')}j (kr)-^j 2 (kr)\, (18) 



where D(z) is the linear growth function and the ji (x) are spherical 
Bessel functions of order I. The first term in brackets accounts for 
galaxy bias while the second describes velocity effects. The ratio 
D(z') I 'D(z) accounts for the growth of perturbations between z' 
and z. The factor dx a /dz converts from Lya flux to the coupling. 
Each resonance contributes a differential coupling (see eq.|6j 



dj, 



(n) 



(19) 



dz' dz' ' 

with the differential comoving flux in Lya from equation 1 141 . 

Because this correlates with the density field, it is easiest to 
observe via 



P M 2(fc) =2P s (k) 



Xtot 



(20) 



The first term probes fluctuations in Tk and ki_o (all encoded in 
p). We show P M 2 in Figure|5| 5 contrasting cases that include only 



5 Our results for J ' P^i are a factor of V2w 2 M 4.4 greater than those 

of BL05, who made an error in their Ps(k) normalisation (R. Barkana, 
private communication). Note that all of their density-induced fluctuation 
amplitudes should increase by a similar amount. 
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Figure 5. Top panel: P 2 power spectrum, which best illustrates the Tj, 
fluctuations arising from density-sourced x a . From bottom to top, the cases 
include only photons from Lyct to Ly/3, from Lya to Ly5, and for n < 
23. Dashed lines indicate / rocyc i c = 1 while solid lines use the proper 
conversion factors. The dotted line is 2f3Pg . Vertical lines indicate the scales 
corresponding to the Lya horizon r a and the HII region size rmi. Bottom 
panel: Transfer function W(k). 
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Figure 6. Ratio of W(k) calculated using proper atomic physics to 
W(k)with / rcC y C i = 1. From left to right, vertical lines indicate the scales 
associated with the Lyct and Ly/3 resonances and with f'Hll- W(k) displays 
small amplitude ripples, which arise from integrating an oscillating kernel 
(the spherical Bessel functions in equation llSI over finite extent (the Lyn 
horizons). Changing /recycle modifies the phase of these ripples leading to 
the wiggles seen in W rccyc \ c (k) /W (k) on intermediate scales. 



photons with v a < v < vg, u a < v < Us, and the entire Ly- 
man continuum. For the latter two models, we show results with 
/recycle = 1 and with the proper atomic physics. Note that each is 
separately normalised to x a — 1. The dotted curve isolates 2P&f3, 
which clearly dominates on small scales. Note that we have applied 
two cutoffs to the power spectrum in this regime (BL05). The first is 
due to baryonic pressure, which prevents collaps e on small scales. 
The se cond is the thermal width of the 21 cm line. lNaoz & Barkanal 
(2005) showed that this thermal cutoff displays a characteristic an- 
gular dependence, which, theoretically, allows fluctuations and the 
cutoff to be separated. To allow easy comparison with BL05, we 
do not include this refinement. Additionally, we expect power from 
the HII regions surrounding the sources to become important on 
scales smaller than the size of a typical HII region thii- We have 
marked this scale for a n isolated galaxy in Figure [5] but note that 
iFurlanetto et alj ( 2004 ) predict that the HII regions could be a fac- 
tor of a few larger at these early times. On sufficently large scales 
kr « 0, the second term in equation 1201 dominates, and W(k) is 
fixed by the source bias. 

Figure|5|cleaiiy shows that the Lyn resonances are important 
on intermediate scales. On large scales only the average flux mat- 
ters, but as we move to smaller scales the higher-n levels become 
important. Figure|6|shows that the fractional reduction in W(k) on 
small scales is ~ 0.63. Although this is near /recycle for Pop. Ill 
stars, that is not the origin of this scaling. In equation 1 1 8i . we can 
write 

%a — ^ &n /recycle (^) : (21) 
n 

and 

W ( k ) = — H & ™ A-yde (n), (22) 



where the b n and a n are defined by reference to equation 1 1 8i and 
the integral of equation I19i respectively. The former care only 
about the local flux, but the latter are averaged over the entire Ly- 
man continuum. If /recycle = constant, they cancel out of W(k) 
and are relevant only as an overall normalisation of x a . In actuality, 
/recycle is a function of 7i, reducing the power. We stress that this is 
because the /recycle (n) are essentially frequency dependent and so 
distort the flux profile about any isolated galaxy (see Fig.|4j- 

It might be hoped that the discontinuities in Figure |4] would 
leave a clear feature on the power spectrum, especially one associ- 
ated with the loss of all photons entering the Ly/3 resonance. Such a 
feature, whose angular scale would be determined by simple atomic 
physics, could set a standard ruler that could be used to test varia- 
tions in fundamental constants or to measure cosmological param- 
eters. Sadly, as can be seen from Figure|5| there is no truly distinct 
feature. Still, the power does decline around k a and measuring its 
shape can constrain the angular diameter distance: we find that the 
amplitude of P i changes by a few percent if the angular diameter 
distance changes by the same amount. However, such constraints 
would also require the astrophysical parameters to be known pre- 
cisely, which will be difficult. 

6.2 Poisson fluctuations 

We turn now to brightness fluctuations uncorrelated with the under- 
lying density perturbations, which can be extracted from the power 
spectrum because of the redshift space distortions (BL05). Specifi- 
cally, if the number density of galaxies is small, then Poisson fluctu- 
ations can be significant. To calculate the correlation function from 
Poisson fluctuations, we again follow BL05 and consider the Lya 
flux from sources within a volume element dV at two points A 
and B separated by a comoving distance I. The correlation function 
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Figure 7. Top panel: Correlation function for the Poisson fluctuations. 
Line conventions are the same as in Figure |5] Vertical lines show 2r a , 
2r@, and 2rnii- Bottom panel: Ratio of the correlation functions assum- 
ing /recycle (") and /recycle = 1- 
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Figure 8. Power spectrum for the Poisson fluctuations. Line conventions 
are the same as in Figure|5] Vertical lines indicate the scales corresponding 
to 2r a and 2rHii- 



takes the form 

xlJv Jm dM r\ r% F ml (z' A )' 

(23) 

where z' B — z '(t\b) is the redshift of a halo at a comoving distance 
tb from a gas element at redshift z. In this expression, we integrate 
over a half volume such that r a < rs , with the factor of 2 account- 
ing for the contribution of sources that are nearer to B. The factors 
P(z') serve to normalise the flux from dV such that 

dx a = P{z')\( M^ff-dMdV, (24) 
v Jm dM 

which makes explicit the expected 1/r 2 dependence of the flux. 
Because of the finite speed of light, points A and B see the sources 
within dV at different stages in their evolution. Following BL05, 
we account for this with the last factor, which scales the source flux 
by the fraction of mass that has collapsed at the observed redshift. 
This ignores a possible dependence of the formation rate on halo 
mass, but in practice high redshift galaxies are highly biased and 
occur within a small mass range just above the minimum cooling 
mass, so this dependence will be weak. Equation I23i is easy to 
understand. For Poisson statistics, the variance of flux from a set of 
identical galaxies would be oc m 2 al n ga il/; this must then simply 
be weighted by the flux reaching each of the two points. Note also 
that, contrary to the claims of BL05, £p oc 1/ fd, where fd is the 
duty cycle of each galaxy, because the fluctuations are weighted 
by two powers of luminosity [MP(z)] but only one factor of the 
density. 

The top panel of Figure^ shows how the correlation function 
increases toward small scales. This is a result of the 1/r 2 depen- 
dence of the flux, which weights the correlations to small scales. 
Including the Lyn resonances amplifies this, because the horizon 
scales skew the flux profile to small radii (see Figure[4}- On large 
scales the correlation function decreases as the two points A and B 
share fewer sources. For r > 2r a a single source cannot affect both 



points so £p = 0. Including /recycle (n) reduces £p, especially on 
the smallest scales, because it decreases the efficiency of coupling 
from level n. The bottom panel of FigureQshows that the suppres- 
sion on small scales is (0.63) 2 « 0.40, because £p depends on two 
powers of the flux. Note that, on large scales, £p increases with the 
proper /recycle- This results from the way we have normalised to 
x a — 1, which reduces the flux of a given source on large scales as 
/recycle increases (see Figure|4j- As with W(k), the scale depen- 
dence is weak except on large scales, where rapid changes occur 
at the appropriate horizons. Proper treatment of the recycling frac- 
tions is clearly necessary to understand the shape of £p. 

These features have similar effects on the power spectrum of 
fluctuations uncorrected with the density fluctuations 

P w _ e (k) = P M o-^-=(^-) a pp(k), (25) 

4 -<> 4 \ x tot J 

as shown in Figure|8| On large scales, taking /recycle = 1 slightly 
amplifies the power. On small scales, they significantly reduce the 
power by ~ 60%. They also affect the shape of the power spec- 
trum, especially near k — n/r a , where the lack of Ly/3 — ► Lya 
imprints a knee on the power spectrum. We also note that the sharp 
Lyn horizons imprint weak oscillations on the power spectrum, es- 
pecially if /recycle = 1 (though these will likely be smoothed by 
photon diffusion). 

6.3 Nonlinearities in the Wouthuysen-Field coupling 

To this point, we have used equation to compute the bright- 
ness temperature fluctuations. This assumes that all the underly- 
ing perturbations are linear; obviously, at x a = 1, this may only 
be marginally satisfied. When the radiation background is large, 
the brightness temperature becomes insensitive to the coupling 
strength and Pr b will be smaller than our estimate. When will such 
corrections become important? One obvious test is whether the typ- 
ical fluctuation Tt is comparable to the maximum brightness tem- 
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perature decrement between coupled and uncoupled gas, 8T (i.e., if 
T s — Tk in eq. 1). But nonlinearities may be important even if this 
condition is not satisfied. A universe with discrete strongly coupled 
regions separated by uncoupled IGM could have small rms varia- 
tions, even though nonlinearities are extremely important in fixing 
the brightness temperature of the strongly-coupled regions. 

Instead we must look deeper at the nature of the fluctuations. 
First, note from Figure [4] that individual galaxies most likely pro- 
vide only weak coupling: x a -C 1 except near to the sources, at 
least if small galaxies (near the atomic cooling threshold) are re- 
sponsible for most of the radiation background. This is not surpris- 
ing: because of Olber's paradox, each logarithmic radius interval 
contributes equally to the background flux in a homogeneous uni- 
verse. The higher Lyman-series photons, together with the finite 
speed of light, also do not dramatically increase the weighting on 
nearby radii. Thus we expect a substantial fraction of the flux to 
come from large distances, where density fluctuations are weak. 
This immediately suggests that the density-dependent power spec- 
trum described in Section loTl will not require substantial nonlinear 
corrections. 

More quantitatively, the fluctuations become nonlinear when 
S Xa = W(k)8(k) > 1; thus we require 

W(k)>l/[a(R)D(z)}, (26) 

where a (R) is the typical density fluctuation on scale R ~ 1/k. 
Figure shows, however, that W(k) is of order unity only for 
k < 0.1 Mpc -1 , where the density fluctuations are themselves 
tiny at these redshifts. Thus, we conclude that a linear treatment 
is adequate for computing the P^i power spectrum, because it is 
primarily driven by large scale fluctuations. 

The Poisson fluctuations in Section l6~2l are more problematic. 
By definition, 

(,p{ta - r B ) = (x a (r A ) x a (r B ))/x 2 a - 1. (27) 

Thus, on scales at which £p > 1, the radiation background near 
galaxy overdensities on this scale is considerably larger than its 
average value, indicating that nonlinear effects are important. In the 
particular model we have examined, £p is large only on comoving 
scales < 10 kpc, so nonlinear effects are again negligible. However, 
the amplitude of the Poisson fluctuations increases rapidly as the 
source density decreases: in models with fewer sources at x a = 
1, or which strongly weight massive galaxies, nonlinearities may 
be important. A maximum value to the variance on any scale is 
ST 2 Q(1 - Q) < 0.25 ST 2 , where Q is the volume filling factor of 
regions with x a > 1. In particular, a linear treatment for sources 
with f-i oc m 2,/3 can violate this limit on scales k > 1 Mpc (see, 
e.g., Figs. 5 and 7 of BL05); in those cases the observed fluctuations 
can be much weaker than linear theory predicts (though it will also 
be non-gaussian). 



7 CONCLUSIONS 

In this paper, we have explored the effects on the spin-kinetic tem- 
perature coupling of photons that redshift into Lyman resonances. 
First, we considered the effect of direct coupling via resonant scat- 
tering of Lyn photons. We showed that the possibility of cascades 
greatly reduces the number of times a Lyn photon scatters before 
escaping. Consequently, the coupling is negligible and may be cor- 
rectly ignored. A side effect of the reduced scattering rate is to make 
the Lyn contribution to IGM heating extremely small, even if T n is 
not in equilibrium with Tk ■ 



Next we considered the increased Lya flux that results from 
atomic cascades. Following the selection rules and transition rates, 
we calculated the probability that a Lyrt photon is converted into 
Lya and showed that /recycle — * 0.36 as n increases. This is sig- 
nificantly smaller than the value /recycle = 1 usually assumed [e.g., 
by BL05]. For a typical Pop. Ill source spectrum, we showed that 
only 63% of the emitted photons will be converted into Lya, de- 
laying the onset of coupling (x a = 1) for a given set of source 
parameters. 

Incorporating the correct /recycle modifies the flux profile of 
an individual source and reduces the coupling on small scales by 
about a factor of 3 for fixed source parameters. In addition, the cas- 
cade process imprints discontinuities onto the flux profile. Using 
the correct atomic physics reduces the amplitude of these disconti- 
nuities and removes one due to the Ly/3 resonance. Unfortunately, 
their weakness is likely to frustrate attempts to use these disconti- 
nuities as a standard ruler. 

We then recalculated the power spectra of BL05, incorporating 
the correct /recycle- This showed a reduction in power of ~ 37% 
(Figure[5J on intermediate scales for density correlated fluctuations 
and of ~ 64% (Figure[S) on small scales for fluctuations uncorre- 
cted with the density. It is possible to mimic this loss of power by 
changing the shape of the stellar spectrum. On small scales, a re- 
duction in the star formation rate will produce a similar reduction in 
flux. Incorporating the proper / reC ycie (n) is thus crucial to correctly 
interpreting 21 cm observations near the time of first light. On the 
other hand, the effects that we have described become unimportant 
once Lya coupling saturates so that T s — > Tk. In this regime, fluc- 
tuations in the Lya flux have little effect on the power spectrum and 
perturbations in the density and reionization fraction dominate. 

Several experiments, including LOFAR, MWA, PAST, and the 
future SKA 6 , are aiming to detect 21cm fluctuations of the type 
we have discussed here. Hopefully, they will be able to study the 
first sources of light through their effect on the IGM around them. 
The details of Lya coupling will determine the observability of this 
epoch. 

We thank M. Kamionkowski for helpful discussions and also 
our anonymous referee for several helpful comments during the re- 
vision process. This work was supported in part by DoE DE-FG03- 
920-ER40701. 
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are terminating hypergeometric series 

2n'P_ = (-l) r [(n - n') 2r {2l > + r)!/(2i>)!] 

x [(n + n')F(-r, -n + /> + 1; 2Z> + 1; v) 

- {n - ri)F{-r,-n + l>;2l> +l;v)], (Al) 

2nP+ = (-iy [(n-n') 2r (2l> +r)!/(2!>)!] 

x [(n + n')(n-l > )F(-r,-n + l > + 1; 2i> + l;v) 
- (n-n')(n + «>)F(-r, -n + Z>;2Z> + l;v)]. (A8) 

These expressions can be inserted into equations {5J and i\2\ to 
compute Pif and /recycle They are in good agreement with ex- 
isting experimental measurements 7 . 



APPENDIX A: CALCULATING EINSTEIN A 
COEFFICIENTS IN THE HYDROGEN ATOM 

The simplicity of the hydrogen atom permits us to compute matrix 
elements for radiative transitions analytically. We are in terested in 
the E instein A coefficients, which may be written as ISobelmat] 
1972) 



A(n,l,j,n',l',f) 
64tt 2 



3h\ 3 



e(2f + 1) 



1/2 
1 



(Al) 



Here {. . .} denotes the Wigner 6 — j symbol and we assume spin- 
1/2 particles. In this expression, the matrix element R™ f has the 
usual quantum numbers n, I, and n r with — 1 = n r . We use the 
sets (n, n r , I) and (n' ; n' r ,l') to describe the upper and lower levels 
in the spontaneous transition. In addition, we let /> be the greater 
of I and V . Because they are separated by Al — I' — I = ±1, 
2/> =1 + 1' + 1. 

We next define, for n' — i> > 1, 



r = ri - l> - 1 > 0. 



(A2) 



We also let 



u — (n — n')/(n + n); 



(A3) 



v = 1 — 1/u 2 = — Ann' /(n — n') 2 



(A4) 



— v/(v — 1) = Ann /(n + n') 2 



-jTi ,1 



(A5) 



Thus we only n eed R„ ', = a pR, where ao is the Bohr radius. 
This is given bv <RudniclJl9'3l 

R = 2 i+i +V + V i+3 {n ~ n) 



(n + n ') n+n ' +1 
(n + l)\ 



-,1/2 



P±. (A6) 



(n' + l')\{n - I -l)\{n> -V 
Here P- and P+ (the subscripts corresponding to the sign of Al) 



7 See reference data at the National Institute of Standards and Technology, 
http://physics.nist.gov/ 
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